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Recent numerical simulations have shown that the distribution of conductances P{g) in three 
dimensional strongly localized systems differs significantly from the expected log normal distri- 
bution. To understand the origin of this difference analytically, we use a generalized Dorokhov- 
Mello-Pereyra-Kumar (DMPK) equation for the joint probability distribution of the transmission 
eigenvalues which includes a phenomenological (disorder and dimensionality dependent) matrix K 
containing certain correlations of the transfer matrices. We first of all examine the assumptions 
made in the derivation of the generalized DMPK equation and find that to a good approximation 
they remain valid in three dimensions (3D). We then evaluate the matrix K numerically for various 
strengths of disorder and various system sizes. In the strong disorder limit we find that K can be 
described by a simple model which, for a cubic system, depends on a single parameter. We use this 
phenomenological model to analytically evaluate the full distribution P{g) for Anderson insulators 
in 3D. The analytic results allow us to develop an intuitive understanding of the entire distribution, 
which differs qualitatively from the log-normal distribution of a Q ID wire. We also show that our 
method could be applicable in the critical regime of the Anderson transition. 

PACS numbers: 73.23.-b, 71.30., 72.10. -d 



I. INTRODUCTION 

The full distribution of conductances P{g) for non- 
interacting electrons at zero temperature has recently 
been studied in detail in quasi one dimension (QID) both 
analytically-'^ and numerically^''^. Large mesoscopic fluc- 
tuations lead to several remarkable features in the dis- 
tribution, including a highly asymmetric 'one-sided' log- 
normal distribution at intermediate disorder bet^ween the 
metallic and insulating limits^, and a singularity in the 
distribution near the dimensionless conductance g I in 
the insulating regime^. While some numerical studies ex- 
ist for 3D finite size systems^iSiSiiSiii, there is no analytic 
method currently available to study the full distribution 
P{g) in 3D, especially at strong disorder. Theoretical 
■work based on 2 -I- e dimensions (e <C 1), 'where a -weak 
disorder approximation can be applied, has been used 
to propose that the critical distribution at the Ander- 
son transition point has a Gaussian center "with po'wer 
la'w tailsiSii^, but this can not be compared 'with numer- 
ical results in 3D2i^ that sho'w a highly non-trivial asym- 
metric distribution similar to the one-sided log-normal 
form of QID. It has not been possible so far to study 
analytically even the simpler case of P{g) in the deeply 
insulating regime in 3D, 'where numerical results point 
to non-trivial deviations from the expected log-normal 
formii. 

The P{g) in QID systems 'were studied analytically 
"within the transfer matrix approach^. In this paper 
'we use a recently proposed generalization of the QID 
approach^^ to obtain analytically for the first time the 
full P{g) for strongly disordered 3D systems. A brief 
account of the "work has been published earlier—. In Sec- 
tions II and III we revie"w briefly the DMPK equation and 



its generalization, respectively. In Sect. IV we analyze 
in detail the numerical data for 3D disordered systems 
in all three transport regimes: metallic, insulating and 
critical. Numerical data allo'w us to determine the free 
parameters of a matrix K which appears in the gener- 
alized DMPK equation. In Sect V we use them to for- 
mulate a simple model for K, and solve the generalized 
DMPK equation analytically. In Sect. VI, an analytical 
formula for the conductance distribution P{g) is derived 
in detail. In our model, the form of P{g) is determined 
by t'wo parameters, T, "which measures the strength of 
the disorder, and 712, "which determines the strength of 
the interaction term in the generalized DMPK equation. 
712 = 1 in the QID systems. The fact that 712 < 1 in 
3D makes the statistics of the conductance in 3D differ- 
ent from that in QID. Although we introduced t'wo ne"w 
disorder dependent parameters, they turn out to be re- 
lated to each other and we sho"w that the present model 
is not in contradiction "with the single parameter scaling 
theory^®. In Sect. VII "we compare the analytical formula 
for the conductance distribution with the numerical data 
and analyze how the distribution P(ln g) depends on the 
parameter 712. In the limit 712 — > 1, we recover the QID 
results. Sect. VIII discusses the possible extension of 
our solution to the critical point. We show that our re- 
sults describe the critical regime qualitatively correctly, 
including the non - analyticity of the critical conductance 
distribution. Finally, summary and conclusions are given 
in Sect. IX. 
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II. THE TRANSFER MATRIX APPROACH 

The distribution of conductances for non-interacting 
electrons at zero temperature can be studied within the 
transfer matrix approach. In this approach, a conductor 
of length Lz and cross-section L x L is placed between 
two perfect leads; the scattering states at the Fermi en- 
ergy then define N o< I? channels. The 2iV x 2A^ transfer 
matrix M relates the flux amplitudes on the right of the 
system to those on the lef1»^. Flux conservation and 
time reversal symmetry (we consider the case of unbro- 
ken time reversal symmetry only) restricts the number of 
independent parameters of M to N{2N + 1) and M can 
be written in general asiiiS 

where m, v are N x N unitary matrices, and A is a di- 
agonal matrix, with positive elements A,;,i = l,2,...iV. 
Microscopic distribution of impurities will lead to a dis- 
tribution pL^{M)d^(M) of the transfer matrices where 
dfj,{M) is an invariant measure which we rewrite as 

PL,{M)dn{M) ^ PLAKu,v)dn{X)d^{u)dn{v). (2) 

If we know the marginal distribution 



PLA{>'a}) ^ / PLA\u,v)dn{u)d^i{v), 
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then the distribution of conductances P{g) can be written 
as 
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Kumar— which has become known as the DMPK equa- 
tion: 
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Here t is the mean free path and the parameter /3 is 
equal to 1, 2 or 4 depending on orthogonal, unitary or 
symplectic symmetry of the transfer matrices. We will 
consider only the case with time-reversal symmetry, for 
which (5=1. Although the parameters Aa are not eigen- 
values of M, it turns out that they determine the eigen- 
values of the matrix TT^ (T is the transmission matrijiSS) 



f*(l + A)- 



(7) 



which characterizes the conductance given by Eq. |(SJ, 
and the matrix v contains the eigenvectors of TT^. So 
we will loosely refer to these as the eigenvalues and the 
eigenvectors in the text. Note that the parameter /3 de- 
termines the strength of 'level repulsion' between eigen- 
values. 

The assumption that pL^{X,u,v) is independent of w, u 
restricts the validity of the DMPK equation to quasi 
one dimension (QID). Quasi ID means not only that 
Lz ^ L where Lz is the direction of the current and L 
is the cross-sectional dimension, but it also requires that 
£, ^ L, where ^ is the localization length. In this limit, 
all channels become 'equivalent', the matrices u and v 
become isotropic and the distribution becomes indepen- 
dent of u or V. The distribution of conductances P{g) for 
such QID systems has been studied recently; it has many 
surprising features arising from large mesoscopic fluctu- 
ations. These include a highly asymmetric 'one-sided' 
log-normal distribution at intermediate disorder between 
the metallic and insulating limits^, and a singularity in 
the distribution near g ^ 1 in the insulating regime^. It 
is not clear if these features persist in higher dimensions. 



is the Landauer conduct ancai^. A systematic approach 
to evaluate the A^-dimensional integral, based on a map- 
ping to a one-dimensional statistical mechanical prob- 
lem, has been developed^, so the full distribution P{g) 
can be obtained if the marginal distribution pL^ii^a}) 
is known. Note that the distribution of other trans- 
port variables which can be written as f{{Xa}), e.g. 
shot noise power^° P = X^^i or conductance of 

N-S (Normal metal-Superconductor) microbridgeSl G = 
Sa=i {i+2\ ' obtained in the same way. 

The above approach is valid in principle for all strengths 
of disorder, in all dimensions. 

If we assume that the distribution p^^ (A, u, v) is inde- 
pendent of M, w, then the evolution of the distribution with 
length Lz can be obtained from a Fokker-Planck equa- 
tion first derived by Dorokhov and by Mello, Pereyra and 



III. GENERALIZED DMPK EQUATION IN 
HIGHER DIMENSIONS 



To study 3D systems, a phenomenological generaliza- 
tion of the DAiPK equation has recently been proposed 
in which the QID restriction is lifted in favor of an un- 
known matrix 
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where the angular bracket represents an ensemble aver- 
age. In terms of this matrix, the marginal distribution 
Pi^(A) satisfies an evolution equation given by^** 
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In QID under the isotropy condition, the matrix K re- 
duces to 



r^QlD _ l±Sab, QID _ -, 

— ' ' lab — 



ab 



N + 1 



(10) 



and one recovers the DMPK equation (with /3 = 1). In 
3D, K is not known analyticaUy, and must be obtained 
from independent numerical studies. 

There are two major assumptions made in Ref. 14] in 
deriving Eq. ©: 

(i) the elements kab can be replaced by their mean val- 
ues Kab, and 

(ii) the L^-dependence of Kab is negligible. 

These assumptions need to be verified before the equation 
can be used. Note that the matrix K depends on the 
choice of representation. Since the assumptions are most 
natural in the position representation, we will study the 
matrix in this representation. 



IV. NUMERICAL DATA 

The generalized DMPK equation apparently intro- 
duces a large number of new parameters, elements of the 
matrix Kab- There is no theoretical prediction about how 
these parameters should depend on the size of the system 
or on disorder. We only know that in the QID limit they 
should follow Eq. (|1U|I . Therefore our first goal is to study 
numerically various 3D and QID systems systematically 
in detail in order to answer the following questions: 

Ql: Are the assumptions (i) and (ii) discussed in Sec- 
tion III valid at all strengths of disorder? 

Q2: How do the elements Kab depend on disorder and 
on the system size? 

Q3: How do the elements Kab depend on the indices a 
and 6? 

Q4: Given the size, disorder and index dependence of 
Kab, is it possible to construct a simple model of 
K at all disorder with only a small number of pa- 
rameters? 

We will address all of the above in this section, but let 
us first briefly discuss the numerical procedure used to 
evaluate K. 

We consider the tight binding Anderson model defined 
by the Hamiltonian 



(11) 



In Eq. 1)11(1 . n = {xyz) counts sites on the simple cubic 
lattice of the size L x L x L^, and e„ are random ener- 
gies, uniformly distributed in the interval [—5,^]. The 
parameter W measures the strength of disorder. The 
Fermi energy is chosen as Ep = 0.01. The hopping term 
tnn' between the nearest-neighbor sites nn' is unity for 
hopping along the z direction and tnn' = t for hopping 
in the x and y directions. Then the dispersion relation is 



E — 2 cos k, -f 2t cos k^. 



2t cos fc„ 



(12) 



For a given cross section of the sample: — Ly — L, 
kx and ky possess values 7r/(L + 1) x 1, 2, ... L (we con- 
sider hard wall boundary conditions). At fixed energy 
E, given values of k^ and ky determine the value of kz, 
which is either real (if jcosfc^j < 1) or imaginary. The 
latter case corresponds to closed channels which do not 
transmit current in perfect leads. To avoid these closed 
channels, which are not considered in the DMPK formu- 
lation, we use t = 0.4. Then the model Eq. ((ll|l exhibits a 
metal-insulator transitions^ at Wc ~ 9. To obtain trans- 
3ort properties, we use the transfer matrix method of 



por t 



[2J|. The main difference from previous works^"-^^ is that 
we also calculate eigenvectors of the matrix TT^ . Using 
Eq. JQ), we calculate numerically the matrix TT^ . Ow- 
ing to Eq. lO, diagonalizing TT^ gives us A as well as all 
elements of the matrix v. 

Note that the eigenvectors depend on the representa- 
tion. In the original formulation of the DMPK approach, 
semi-infinite leads consist of mutually independent and 
equivalent ID wires. Therefore, the transfer matrix in 
the leads is diagonal in both channel and space repre- 
sentations. In numerical work, we need to distinguish 
between these two, since the transfer matrix is diagonal 
only in the channel representation. We therefore calcu- 
late the matrix TT^ in the channel representation, find 
eigenvalues and eigenvectors, and transform the latter 
back to the space representation to obtain the matrix v. 
Elements of v are then used for the calculation of the 
matrix K in the space representation. 

We now go back and address the questions raised at 
the beginning of this section. 



A. Ql: Validity of the assumptions (i) and (ii) 

In order to check assumption (i), i.e. if the elements 
kab can be replaced by their average Kab in 3D, we an- 
alyze the probability distribution P{kab)- We start with 
the weakly disordered metallic regime. Two values of dis- 
order were used: W = 2 and W = A. More relevant than 
the actual strength of the disorder is the mean free path 
£ which can be estimated from the mean conductance^^ 



(g) 



N 



-I 



(13) 



with N = . From the dependence of {g) in QID 
systems we estimate 1{W = 2) « 9.2 and 1{W = 4) w 
1.8, in units of the lattice spacing. 
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FIG. 1: Probability distribution of normalized kab for QID 
systems, (o: a = fe = 1, 0: a = 1, 6 = 2, V- a = b = 2). 
Full symbols: W = 2, L = IQ, Lz = 16L; open symbols: 
W = 2, L ^ 10, Lz = 24L. Lines with symbols: W = 2, 
L = 14, Lz — SL. Data confirm that the distribution becomes 
narrower when L increases. Inset shows the same for W = 4 
and L = 18, = 4L. As expected, distributions are broader. 
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FIG. 2: The distribution of normalized kab for 3D systems 
with ly = 4 and various system sizes, (critical disorder Wc ~ 
9.) Open symbols: a = b = 1; full symbols: a = 1, 6 = 2. 
Inset shows the L-dependence of Aab = \/var kab/ Kab which 
decreases when L increases. Data confirm that distribution 
is self- averaging (it becomes narrower and Aab — > when L 
increases) although it is broader than in the QID case. 
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FIG. 3: Probability distribution of normalized matrix ele- 
ments kab- (a) P{kaa/Kaa) = 1,2 aud 3) at the critical 
point W — Wc ~ 9. Distribution is L-independent (apart 
from the exponential tail which is broader for larger L since 
mean value Kaa ~ ^/L). Note the logarithmic scale on the y 
axis, (b) Distribution of off-diagonal elements ki2 and ^23 
possess sharp maxima close to zero, and long exponential 
tails. Insets show standard deviations of distributions as a 
function of mean values for 8 < L < 18. These data also 
show the accuracy of our estimate of the critical point since 
we expect both the mean and the standard deviation to be L- 
independent at W = Wc- The distributions for QID systems 
L X L X 8L are almost identical to those for cubes (data not 
shown) . 



First, we test how assumption (i) is fulfilled in the 
weakly disordered limit, where the DMPK equation is 
known to describe the universal features of transport 
statistics correctly. We show in Figure ^ the distribu- 
tion P{kab) in the weakly disordered QID regime. As 
expected, the width of the distribution increases when 
W increases, but P{kab) is self-averaging (it becomes 
narrower when L increases). Figure |21 shows the same 
distribution for 3D systems. The distribution is again 



self-averaging, although much broader than in the QID 
case. 

In the critical regime, the probability distribution 
P{kab) is no longer self-averaging but tends to be L- 
independent in the limit L ^ 00 (Fig. |3J. Although 
the distributions possess long exponential tails, they have 
well defined sharp maxima, which do not depend on the 
system size. 

In the insulating regime, the distribution P(fcii) be- 
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FIG. 4: Probability distribution P{kii) in the insulating 
regime {W — 16). Although distribution becomes narrower 
when system size L increases, it is not self-averaging. Inset 
shows the size dependence of the An — Vvar kn/Kn for 
W = 16 and W = 20. An converges to a non-zero constant 
when L oo. 
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FIG. 5: At the critical point, which depends on the 
anisotropy parameter t, both LKi\ and LK\2 converge to 
I/z-independent values when L^jL oo. The limiting value 
LKii depends on t. Shown are also data for the insulating 
regime W = 20, t = 0.4. All data for K12 almost coincide so 
that they are not distinguishable in the figure. 



comes narrower when L increases (Fig. 0}. However, on 
the basis of our numerical data we conclude that the 
distribution is not self-averaging. Although var fcn de- 
creases when L — > 00 (data not shown), the normalized 
width All = Vvar kn/Kn (shown in inset of fig. ^ 
slightly increases when L increases. As Kn itself is non- 
zero in the limit L 00 (fig. EJ, An should converge to 
an L-independent function for large L. 

The distribution of off-diagonal elements P(fci2) in the 



insulating regime (not shown) is qualitatively the same 
as that at the critical point. 

We conclude that both in the critical and localized 
regimes the distributions converge to L independent func- 
tions with a well-defined peak, but the standard deviation 
is of the same order of magnitude as the mean. We note 
that the most-probable value of fcn is always very close to 
its mean value; we therefore expect that replacing fcn by 
its mean value Kn is a reasonable approximation as long 
as one is interested in qualitative results only. Thus, al- 
though to leading order assumption (i) remains valid for 
all disorder, we have to keep in mind that fluctuations of 
the elements kab in the strongly disordered regime might 
become important if the final results are sensitive to the 
exact values of these elements. We have checked that 
the final distribution of conductances do not change in 
any appreciable way if fluctuations of fcn are included by 
overaging the conductance distribution over P(fcii). 

To check assumption (ii), namely if the Lz dependence 
of Kab is negligible, we studied the L^/L dependence of 
Kii and K12. Figure [S] confirms that for both critical 
and insulating regimes the parameters Kn and if 12 con- 
verge to non-zero (although t - dependent) limits when 
the length of the system increases. It shows that the 
properties of the matrix K depend only slightly on the ra- 
tio Lz/ L and reach Lz-independent limiting values when 
Lz/L 00 'm all transport regimes. The assumption 
(ii) is therefore reasonably well satisfied at all disorder as 
long as Lz > L. 

Thus we conclude from our numerical studies that to 
leading approximation the generalized DMPK equation 
(|5J) remains qualitatively valid at all disorder in 3D, but 
the effect of fluctuations of kab on the final results has 
to be evaluated in more detail before a quantitative com- 
parison with numerical results can be made. 



B. Q2: Disorder and size dependence of Kab 

We start with the weak disorder regime. To distinguish 
the generic W dependence of Kn from finite size effects, 
we analyzed in Fig.^lthe L dependence of the parameter 



n= {N + l)Kn/2. 



(14) 



As expected, k decreases when L increases. However, 
from the analysis of QID systems (also shown in Fig.EJ 
we conclude that n converges to 1 only for very small 
values of disorder. As shown in Fig. El ~ 1.36 for 
W = A. Thus, deviations from Eq. H1U|) already appear 
in the metallic limit, probably due to the decrease of the 
mean free path. 

In the next step we analyze how the matrix K changes 
when the disorder increases. Typical results are given in 
Fig. Our data show that the L dependence of the pa- 
rameters LKii and 712 is different for different transport 
regimes. While in the metallic regime LKu decreases as 
~ 1/L and 712 converges to unity when L increases, qual- 
itatively different behavior is obtained at strong disorder. 
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FIG. 6: dependence oi k = {N + l)i5:ii/2 (open sym- 

bols) and (A'' + 1)7^12 (full symbols) in the metallic regime. 
Data confirms that k depends on L. This agrees with data in 
Fig. 1121 Also, 27^12 differs from A'n for small L. This agrees 
with data in Fig. 1131 To estimate limiting behavior of Kn, 
we also considered QID systems L x L x 4L. k is close to 1 
only for very weak disorder W — 2. For Vl^ = 4 we obtained 
K ~ 1.36. As this value does not depend on L for 8 < L < 18 
(0), we expect that 3D data for W = 4 will converge to the 
same value when L <x. 



In the insulating regime Ku converges to a non-zero L- 
independent constant when L — s- oo (Fig. |5J and K12 
converges to zero as K12 ~ 1/L (Fig. I1(J|) . This means 
that 712 ~ 1/L in the insulating regime (Fig. |SJ). 

Figures 7 a,b also show that there exists a critical dis- 
order W — Wc where both LKn and 712 are independent 
of L. Note that 712 < 1 at 1^ = Wc- We found that the 
critical value 712c depends on the anisotropy (Fig.|SJ). For 
the present case t = 0.4, 712c ~ 0.28. The qualitative L 
dependence in different transport regimes is summarized 
in Table n 

We observe that the disorder dependence of Ku is con- 
sistent with Kii oc where m = 2, 1, in the metal- 
lic, critical and insulating limits, respectively, in agree- 
ment with Ref. 26]. Note that in contrast, Ku oc 
for all strengths of disorder in QID. This is a major dif- 
ference between QID and 3D. One can understand qual- 
itatively how the L dependence of Ku changes in the 
weak and strong disorder limits on general grounds. If 
all channels are equivalent, we expect the column ma- 
trix Via ^ which satisfies the unitary condition 
Ea=i l^^iaP = 1- This leads to Ku ~ 1/N = l/L^ in 
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TABLE I: Typical L-dependence of the parameters Ku and 
712 in the metallic, critical and localized regimes. 
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FIG. 7: Disorder dependence of Ku and 712 for various sys- 
tem sizes. Note the common crossing point at W = Wc ~ 9. 
Note also that 712 — > 1 for W < Wc and L — > 00, as ex- 
pected from DMPK, but 712 decreases with the system size 
for W >Wc. 



the metallic limit. On the other hand, if the localization 
length ^ ~ 1, then on any cross-section at a given L^, 
we expect only a few sites on the back side of the sample 
to be 'illuminated' by an incoming wave, so we expect 
't'la ^la- This leads to Ku 1, independent of L. 
Similarly, since all Kia cx in the metallic regime, 

we expect 712 '--^ 1 in the metallic regime. However, in 
the insulating regime, we have not found a simple phys- 
ical argument why K12 (xl/L and hence 712 oc 1/L. We 
also find numerically in the insulating regime that for 
1<C^^L, The structure of the eigenvec- 

tor Via that gives rise to Ku ^ 1/C and K12 ~ 1/L in 
the region 1 <C ^ ^ L is highly non-trivial, and deserves 
further study. 

Figure 1^1 confirms our claim that Ku > in the insu- 
lating regime. Its limiting value, Ku — lim^^oo Kii{L) 
can be used as an order parameter for the scaling analysis 
of the Anderson transition. It is evident that Ku = for 
W <Wc and Ku > for > Wc- We show in Fig.UIl 
the L-dependence of LKu and of the mean conductance 
(g) for various strengths of disorder. Similar behavior is 
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FIG. 8: 712 as a function of the system size for various 
strengths of disorder. In the metallic regime, 712 converges to 
1 when L — > 00. At the critical point, 712 is L-independent, 
and in the insulating regime 712 decreases when L increases. 
Solid lines are linear fits 712 = a + b/L with a ~ 10~^ for 
W < 40. 
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FIG. 10: System size dependence of 7^12 for various values of 
disorder in the insulating regime. Data confirm that K12 ~ 
1/L and that limi^oo K12 = 0. 




FIG. 9; System size dependence of Ku for various values of 
disorder. Solid lines are linear fits. Data in legend presents 
limiting values Ku = lim_L^oo Js'ii (i). Shown are also data 
for the metafile regime {W = 4) for which Ku ~ (Table 
I). While Kii = in metal and at the critical point (W — 9), 
it is non-zero in the insulating regime. Thus, it may be used 
as the order parameter of the Anderson transition. 



shown in Fig.|Slfor the parameter 712. One sees that all 
three parameters, (g), LKn and 712, could be used for 
the estimation of the critical disorder Wc^ at which none 
of them depends on the system size. 



C. Q3: Index dependence of Kab 

Again we begin with weak disorder. To compare 3D 
and QID systems, we show in Fig.^lthe parameters Kaa 
as a function of the index a. It is clear that the 3D data 



differ considerably from the DMPK value 2/{N +1). In 
Fig- El we show the ratio 71a = 2Kia/Kii for various L 
and compare it with QID numerical data. It is evident 
that 71 a converges to 1 when the system size increases, 
in spite of the fact that both Kia and Ku differ from the 
DMPK values of Eq. El The convergence is much slower 
in 3D than in QID systems. Also, 71^ converges slower 
for larger a. We conclude that although 3D metals are 
qualitatively similar to QID metals, there are quantita- 
tive differences that need to be explored further. 

In the critical regime. Figures El and El show that 
the a and b dependence of the matrix elements Kab 
can be described by simple functions: Kaa ~ Kn/a^/'^^ 
7ia ''^ \/La}/'^. Although we did not analyze all matrix 
elements in detail, we believe that the data presented 
here support our expectation that all matrix elements 
Kab can be expressed in terms of Ku, K12 and some 
simple function of the indices a and b. 

In the insulating regime we find a simple a-dependence 
of the difference 

Kaa - i^ll « -CV^ (15) 

with |c| ~ 0.04 (Fig. El- We also found, (Fig. El, that 
for a < N/4:, 71^ ~ 1/y^L, very close to the value ob- 
tained at the critical poinljSi. 

There is an interesting correspondence between the a- 
dependence of Kaa and the a-dependence of the param- 
eters Xa defined through the parameters^ 

Xa = sinh^ Xa (16) 

which is summarized in Table HI! 

We find that the index dependence of Kaa can be ig- 
nored in the metallic regime. The a-dependence is more 
pronounced at the critical point and in the localized 
regime. On the other hand, higher channels (a 3> 1) do 
not contribute to the transport either at the critical point 
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FIG. 11: Estimation of the critical point from the mean con- 
ductance (g) and from LKn. Both parameters (g) and LKn 
are L-independent at = Wc- LKn cx L (oc 1/L) in the 
insulating (metallic) regime, respectively. From data we con- 
clude that 9 < Wc < 9.5 and use Wc = 9 throughout the 
paper. 



FIG. 12: a-dependence of Kaa x {N + l)/2 for 3D and QID 
systems with W — 2, (mean free path I « 9.2) and for W = 4 
{I ~ 1.8). Kaa are larger than is predicted by the DMPK 
theory. Good agreement with DMPK is observed only for 
QID systems with very small disorder {W — 2). 
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TABLE II: Index dependence of Kaa obtained in the present 



work compared to Xd 
regimes. 



I in the metallic, critical and insulating 



FIG. 13: a-dependence of 71a for the metallic regime {W — 4). 
Finite size effects are clearly visible so that 712 — 1 only in 
the limit of large system size. As expected, convergence is 
better for QID systems than for cubes. Nevertheless, data 
confirm that 71a converge to 1 in the limit L — > 00 even when 
Kii and K12 do not converge to values assumed by DMPK 
given in Eq. pOjl . 



or in the localized regime, so the actual values of Kat for 
large a and b are not important. Therefore, we conclude 
that the weak index dependence of the matrix elements 
Kab is less relevant for transport properties compared to 
the dependence of the matrix elements on disorder in the 
L —t 00 limit. 

It is also worth mentioning that since we are interested 
only in L-independent quantities at the critical point, the 
a-dependence of any parameter is relevant only for a < L. 
When a becomes comparable to L, we can not distinguish 
the true a-dependence from finite size effects. 



D. Q4: Simple model for K 

Finally we ask the question if it is possible to construct 
a simple model of Kab with only a small number of inde- 
pendent parameters. We just concluded in the previous 
section that the weak index dependence of the matrix at 
strong disorder is not very important. We expect that 
the crude approximations 

Kaa ~ Kn and -fab ~ 712 (17) 

capture the major qualitative features of the matrix K. 

Approximation (|17(l introduces two new parameters, 
Kii and 712. We show that both parameters Kn and 
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FIG. 14: Critical point: We plot La^^^Kaa to show that 
LKaa oc a"^''^ are L-independent for all a < 20 




FIG. 15: Critical point. From the linear dependence 
(a'''Vi)7ia oc a/N we conclude that 71a oc l/a^^^L. {N — 



7i2 are unambiguous functions of the localization length. 
To do so, we calculated the hmiting values of Ku (Fig.|^ 
and of the parameter xi (Fig. I18|l and plot Ku versus 
£, (Fig. I19|l . We observe that the ^ dependence of Ku 
depends also on the anisotropy parameter t (data not 
shown) . 

In the same way, we analyzed parameter 712. As was 
shown in Fig. |S1 we expect 712 ~ in the localized 
regime. Data in Figures 1^ and ITUl support this assump- 
tion. To estimate the disorder dependence of 712 , we plot 
in Fig.|2Dlthe quantity 



I 




FIG. 16: Insulating regime: a-dependence of Kaa for 3D sys- 
tems with W — 30. To see the a-dependence more clearly, 
we plot the difference [Kaa — A'njy^ which behaves as ~ —a 
for small a. Solid line is a linear fit, from which we have 
Kaa = Kii + 0.05 - 0.04/Va for L = 18 and a < 18 
{Kii = 0.44). Right inset presents data for i = 14 and 
various strength of the disorder. Left inset shows Kaa for 
L = 18. 




0,2 0,3 
a/N 



FIG. 17: Insulating regime {W — 30): a-dependence of 71a 
(L = 6, 10 and 14). Similar to the critical point, data shows 
that a^^'^/L x 712 ~ a/N so that 71a ~ 1/(^L). Inset shows 
data for L — 14 and various strength of the disorder. 



hm L712 = 2L —- (18) 

L^oo hniL^oo Ku 



V. SOLUTION OF THE GENERALIZED DMPK 
EQUATION IN THE STRONG DISORDER LIMIT 



which indeed increases linearly with localization length 
in the strongly localized regime. 



While modeling the full jab at the critical point needs 
more careful numerical studies, the insulating limit is 
simpler and provides a test case for the generalized 
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FIG. 18: Insulating regime: Estimation of the localized length 
^(VF) from the linear L dependence x\(V) = const + 1//^. 
(Ai = sinh^a;i). Values of ^ are given in the legend, g ~ 



cosh xi Ri exp— 2L/^. 
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FIG. 20: 1/712 defined by Eq. HISII as a function of ^. Data 
confirms that there is an unambiguous ^-dependence of 712. 
This is important for the formulation of the single parameter 
scaling theory. Dashed line is the linear dependence 712 = 
^/2L, considered in [T5| . 
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FIG. 19: Figure (a) shows how LKn depends on xi for var- 
ious disorder in the localized regime. Data confirm linear 
relation LKn oc xi. Figure (b) shows how limiting values 
of Kii, obtained from the L-dependence of Kii{L) (Fig. |3Jl, 
depend on values of 1/^ fFig. I18( . Data confirms that there 
is an unambiguous ^-dependence of Kn. 
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FIG. 21: Ratio r/(7i2^) = l/2L_K'i2, given by Eq. HHfor a 
cubic system, as a function of disorder for 3D systems in the 
insulating regime W > Wc- 



DMPK equation (8). It predicts that the logarithmic 
interaction between the transmission eigenvalues Aq van- 
ishes as 1/L in the insulating Hmit. In the rest of the 
paper we will test this prediction by evaluating the full 
distribution of conductances in the insulating limit for 
a 3D conductor as described by Eq. (8), using the sim- 
ple approximate model for K suggested by our numerical 
studies, namely 



Kii and w 712 ~ const. — . (19) 



It is useful to introduce another parameter. 



and the ratio 



r 

712 



11 

L £ 
Tl 2LKi2 ' 



(20) 



(21) 
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r l/Lj <C 1 in the insulating regime, both in 3D 
(L — Lz) and QID (L^ ^ L) systems. It measures the 
strength of the disorder. The ratio r/712 reduces to F in 
the QID hmit when 712 — 1. For 3D {L^ — L) we see 
in Fig. 12] that F/712 ~ L£/Lz varies smoothly between 
about 2£ and 3.5^ when disorder increases from Wc to 
infinity. Thus, the strength of disorder is characterized 
predominantly by the parameter F. For definiteness, we 
will use F/712 = 2 appropriate for a cubic system when 
needed for comparison with numerical results^S. 

Note that fluctuations in fcaa, ignored in the model, 
would lead to fluctuations of F and 712 as well, but will 
not change either the length or the disorder dependence 
of these parameters. 

A brief description of the results has appeared in . 
Here we provide many of the details. 

We rewrite the generalized DMPK asi^ 



d{L,/e) 



d 

' dXa 



dp 

dXa 



on 

^ dXa 



(22) 



where 



N N 

J7 = — In I sinh^ x^ — sinh^ Xa^'-'' — In sinh 2xa- 

a<b a—1 

(23) 

We define P = e""/^*. Then following 

'3/ satisfies the imaginary time Schrodinger equation 

-d^/d{Lz/e) = {n- [/)* with 



n 



a<b 



1 

dxl sinh^ 2xa 
1 



sinh^(xQ — Xb) sinh^(a;a + Xb) 



(24) 



and U constant, where the strength of the interaction is 
given by 



V = —;— 712(712 - 2). 



(25) 



The interaction term in QID vanishes for the unitary 
case, when j'^-l^ =[3 = 2. Note that the interaction also 
vanishes in the limit 712 — > 0. In the 3D insulating case, 
712 ^ <C 1, and the interaction can be considered 

negligible, for all symmetries. We can therefore use the 
712 = 2 solution of 30] for our 3D insulators. The solu- 
tion in the insulating limit is then given by'^" P = e^^ , 
with 



N 



a>b 
N r 

-E 



In I sinh^ Xb — sinh^ XaP^^ + In |a;^ — x1\ 



a=l 



In sinh 2xa + lnXa — Tx^ 



The replacement of /3 = 2 in Ref. jS^I by 712 



(26) 
in 



Eq. H26|) has the consequence that while all (xa) 3> 1 in 



the insulating regime, the difference s = (Xa+i—Xa) is not 
of the same order as (xa)- For example, if we keep only 
the first two levels, the saddle-point solutions for xi and 
X2 give (xi) ~ Lz/S, and {x2 — xi) <C (xi). We therefore 
do not assume that X2 3> xi. However, we do make the 
simplifying approximation that In | sinh^ Xa — sinh^ Xb\ ~ 
Insinh^ Xa and In \x'^ ~ ^fl ~ lii^;^ for a > 2 and a > b. 
Eq. H26|l then becomes 

JV 

H^H,+Y, [V{xa) - 7i2(a " 2)/(x„)] , (27) 



a=2 



where 



Hi = - In 1x2 -xi\ + Txf - -lnsinh2xi -Inxi, (28) 



V{x) = Fx^ — — In sinh 2x — In x — 712 In sinh x, (29) 



/(x) = In sinh X + fclnx; k = . (30) 

712 

We can now use the method developed in Refs. to 
obtain the full distribution P{g). 



VI. P{g) IN 3D IN THE INSULATING LIMIT 

As in PQ, we separate out the lowest level xi and 
treat the rest as a continuum beginning at a point X2. 
Then 



H = Hi + I dx'cr(x') 



where we have used 



V{x') - 7i2/(x') / dx"a{x") 

(31) 



(a - 2) = / dx'a{x'). 
The density satisfies the normalization condition 



a{x)dx = N -1. 
P{g) can be obtained from H{{xa}) as 

P{9) = J-JH dxae-"5{g - ^ sech^x,,) 



(32) 



(33) 



(34) 



where the 5-function represents the Landauer formula for 
conductance. It turns out that because of the nonlin- 
ear dependence of the Hamiltonian on the d ensit y, the 
complex delta function representation used in |ll4j is not 
suitable for the present case. We therefore use a real 
representation 



S{x — a) = lim 



1 



(35) 
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Following fl"^, P{g) may be expressed as 



oo poo 



P[g) = / D[a{x)] / dx / da:2e'-^(^^'^^''"(^» (36) 
where the free energy functional F is given by 



F= Hi+ I dx'a[x') 



1 

'A2 



(9 -hi) -I 



where 



and 



V{x') - 7i2/(:e') / dx"aix") 

J X2 



(37) 

j^(^'M.i^'\ (38) 



K = 2 / da;'cr(a;')/i(a;') 



= sech^x; hi — sech^xi. 



(39) 



The saddle point density is to be obtained by minimiz- 
ing F with respect to cr(x), subject to the normalization 
condition. We therefore define 



T = F - K 



a{x)dx -{N -I) 



(40) 



and minimize JF. It is useful to rewrite the free energy 
using the normalization condition in a way that removes 
the upper limit from the resulting equation. We use 



Sa 



dx'a{x')f{x') / dx"cr{x") 

b 



(TV - 1) - / dx"a{x") 



— I dx'a{x')f{x') 



= fix) / dx'a{x') - / dx'a{x')f{x') (41) 
to obtain 

fX px 

Y{x) - A = 712 f[x) / dx'a{x') - / dx' a{x')f{x') 

J X2 J X2 

(42) 

where we have defined 

Y{x)^V{x)-^[{g-h,)-^]h{x). (43) 

Eq. (|42l) evaluated at a; = a;2 fixes A = Y{x2)- Taking a 
derivative of Eq. (|42|l with respect to x (represented by 
a prime) gives 



Y'{x) = 7i2/'(x) / aiy)dy 



(44) 



Evaluated at a; = a;2 , this fixes X2 as the beginning of the 
continuum 



Y'{x2) = 0. 



(45) 



Taking another derivative with respect to x, it is now 
possible to obtain the density 



a{x) = — ' 



712 V /'(^) 



(46) 



We check that plugging a{x) in Eq. back to Eqs. (03 
I44|l satisfy those equations. The density has the form 

1 K 
a{x) ^ ai{x) ~ —a2{x) + ^b{x). (47) 

Plugging this form in the definition of k, wc obtain 

2 h{x)[ai{x) - ■^a2{x)]dx 



1 2 



■ h{x)h{x)dx 
which, expanded in powers of A, is given by 
^t = 2((7-/ii) + AVi+0(A^). 



where 



and 



Ml = ^[-ai + 2(g-/ii)] 



(48) 



(49) 



(50) 



Oil 

(3 = 



712 Jx2 \ / W 



V'{x)\ 



112 J X2 



h'{x)\' 



h{x) 
h{x). 



(51) 



Using the expansion for k, given by Eq. (48), we obtain 
in the limit A from Eq. (42) 

Y{x) ^V{x)+ ^Jiih{x). (52) 

The free energy can then be written as 



F{xuX2) = H,+ f dxa{x)\vix)^ f{x)^j^ 
Jx2 I f [x) 



(53) 



There are two additional constraints that were not in- 
cluded in the variational scheme and will be enforced 
directly. 



a{x2) > 0; X2 > xi 



The conductance distribution now becomes 



(54) 



P{g) = f dxi f da;2e-^("i'"^'S)5(y'(a;2)), (55) 

Jo Jxi 

with equations (j^Hl, lEOl, 63 and defining V{x), 
f{x), a(x), and Y{x), respectively. 
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A. The free energy 



B. The constraints 



Let us write F{xi, X2) ~ Hi {xi, X2)+F2{x2) and define 
W{x) ^Y'{x)/f'{x). (56) 
Then a{x) = W'{x)/ji2 and 

F2^— I dxW'{x)V{x) - I dxW'{x)f{x)W{x) . 

(57) 

On the other hand, defining 

$(x) = V{x) ~ ](x)W{x) (58) 
and using partial integration, we get 
1 

F2 ^ / dxW{x) [V\x) - f'{x)W{x) - f{x)W'{x)] , 

Jx2 

(59) 

where we have neglected an irrelevant term ^{b)W{b)/ J12 
independent of X2 and we have used W{x2) = 0. Using 
f{x)W{x) = Y'{x) and V'{x) - Y'{x) = -/xi/i'(x), we 
rewrite the above as 



We already have one constraint Y'{x2) = 0. We also 
demand cr(x2) > which requires 

Y"{x2) > 0. (66) 

This defines X2min- Also, from Eq. (|5n|l . 

2 1 

g -1/ cosh xi = -(«!+ /ii/3) = go, (67) 



or 



xi = cosh 



1 



(68) 



On the other hand defining Wi{x) = V'{x)/f'{x) and 
1^2(2;) = h'{x)/ f'{x), it is easy to see that 



1 



F2 = — / dxh'{x)W{x) + / dxf{x)W'{x)W{x). 

Jx2 712 Jx2 

(60) 

The two alternate expressions for F2 can now be com- 
bined to obtain 



£dxf{x)W'{x)W{x) 



1 



^ / dxh'{x)W{x) + C, 



dxW'{x)V{x) 



(61) 

where C is a constant. Plugging this back to Eq. ()60|l . 
we obtain 



Fo - 



1 



271 



dxW'{x)V{x) + 



Ml 

2712 



dx/i'(a:)VF(x). 



(62) 

We can again use partial integration to rewrite the first 
term as an integral over V' {x)W{x), using again the fact 
that W{x2) = 0. Then the two terms can be combined 
to obtain 



F2 



We define 



1 



dx 



2712 Jx2 n^) 



^[T/-(x)-M?/.-(.)] 



(63) 



■h ^ dx——^; J2 



J?, — I dx 



fix) ' 
h'^{x) 



" V'{x)h'{x) 

X2 /'(^) ' 



(64) 



Then 



Fo = 



271 



-[Ji-^iiJs]; Ail = -77 



V'{X2) 



h'{x2) 



(65) 



go — dxh{x)a{x). 



(69) 



Using partial integration and the fact that ^'(2:2) = 0, 
we can also rewrite 

go^-—[.h + ^il•h]■ (70) 
712 



C. Validity of the approximations 

We started with the assumption that while (a;i) ^ 1 
in the insulating regime, the difference {x2 — xi) ^ {xi). 
It is therefore important to estimate the difference from 
the above results. We will use saddle points of the free 
energy 

F{xi,X2)f=i V{xi) -\n{xl- xl) + F2{X2)] ^ 

V{xi)^Txl-xi-\nxi ^ > 

where F2 is given in Eq. I|63|) . Let us define s — X2 — xi 
and C — xi ~ 1/2T. We assume s xi and Q xi. 
The saddle point solutions for s and C are obtained from 
dF/dQ — and dF/ds — 0. Using chain rule to write 
the partial derivatives in terms of dF/dxi and dF/dx2 
we obtain 



v'ixiy 



X2 + Xi 



-+F^{x2) = 0; 



1 1 



S X2 + Xl 



-+F^{X2) = 
(72) 



where prime denotes derivatives with respect to the ar- 
guments. Combining the two gives 



From the definition of F2 we have 



(73) 



F^{X2) - - 



2712 



dJi 
dx2 



A*i 



,dJ3 
dx2 



^^Js. (74) 
712 0x2 
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The integrals Ji and J3 depend on X2 only via the lower 
limits of the integrals. Derivatives w.r.t. X2 are simply 
the negatives of the integrands evaluated at the lower 
limit. Using definition of /ii this gives the first term in 
Eq. (O equal to zero, leaving ^^(0:2) = J3. Using 

h{x) « 4e^^^ and 1/ f'{x) Ri x/k for x/k ^ 1, we get 



- [16x26- 



4X2 



4e 



-4X2 



(75) 



Then 



J(i + r2).x2 + ^(i + ri), (76) 



where T3 = r(l - 3r/4), T2 = -QT + and 
Fi = — r/3. We neglect l/(x2 + xi) compared to the 
other terms in F'{s) = (Eq. ifT^ ) and expand ^2(2^2) 
in Taylor series around X2 = l/2r. The dominant term 
is i(C + s)2r where we have used F^"{x2 = l/2r) = T. 
Using ( — — l/2rs, we finally obtain 



As we will show, the actual density of the levels play a 
minor role in the distribution P{g), which is dominated 
by the first few levels. Therefore our results will be qual- 
itatively correct, although there would be quantitative 
discrepancies due to the difference in the density. 

Note that in the opposite limit X2 ^ k, x ^ X2, the 
density becomes 



aix) » 2r/7i 



(81) 



This corresponds to a uniform average spacing s = 
{xa+i — Xa) of eigenvalues of order unity (L = L^), com- 
pared to the uniform spacing s ~ L^/C in QID. In con- 
trast, 3D metals are similar to QID metals having uni- 
form (j{x) extending down to x = and s ^ L^/L'^. The 
opening of a gap in the spectrum of Lyapunov exponents 
~ 1/^ may be considered as the signature 
of the Anderson transition. 



VII. RESULTS AND DISCUSSIONS 



dF 



1 



1 1 



-Es^ ^ 0. 



ds 8Es2 2 ' 2' 
Therefore, the saddle point solutions are given by 



(77) 



(78) 



Since (xi) ^ 1/2E, we confirm our expectation that (x2 — 
xi) <C (2:1 )• The results are also consistent with our 
assumption that both s and C are much smaller than xi, 
so the free energy calculations remain valid. 

However, numerically we find that while (xi) ~ 1/2E, 
s 1 independent of disorder. Thus our result s — 
^ (xi) is only qualitatively correct. The fact that 
actual s is much smaller than what we find is related to 
the inaccuracy in our evaluation of the density. Indeed, 
we can obtain the density directly from Eq. (|46|) . In the 
limit X ^ 1, 712 <C 1 we find 



2F 



+2V'{X2) 



fc(2rfc+i) 

(x+k)^ 



2{x+k) 



k ^ \ ^-2{x-X2) 



(79) 



where we have used k = 2/712 (Eq pOfl l. In the limit 
X2 -^k but X > X2, the density simplifies to 



cr(x) « 2rx. 



(80) 



The linear x-dependence as well as the F dependence 
agrees with numerical results. However, the slope turns 
out to be too large. This is possibly the consequence of 
our simplification of the Hamiltonian Eq. H26|l to Eq. (|27f) . 
where all the interaction terms were neglected except for 
the one between the first and the second levels. As shown 
in Q , it should be possible to obtain an integral equation 
for the saddle point density which can then be solved at 
least approximately. 



With the above caveat in mind, the saddle point free 
energy Fsp{xi,X2; g) has the form (Eq. lO) 



Fsp{xi,X2) = Hi 



1 r dx 



2712 Jx2 fix) 



(82) 

where primes denote x-derivatives. Eq H55|l can then be 
rewritten as 

P{\ng)ocg dxae-^^'^^i'^'^^^^'^e-^C^^-xi)^ (gg^ 



where the integration over xi is eliminated by a con- 
straint arising from the minimization of the free energy: 



ci = cosh ^[l/y/g - go]; 



90 = -^1x2 J^h'{x)[V'ix)+^,lh'ix)]. 



(84) 



The lower limit X2mm is the larger of the additional 
constraints imposed by the conditions cr(x2) > and 
X2 > xi > 0, xi real. 



A. Analytical model 

It is instructive to consider first a simple approxi- 
mate solution of Eq. 1)83(1 . which is dominated by the 
lower limit of the integral. To a good approximation, 
go is negligible compared to g in the insulating limit, 
and xi « iln(4/g). The condition (t(x2) > or 
equivalently Y"{x) > gives X2mm from the condition 
Y"{x2rmn) = 0. Thls gives 



X2r, 



(i + r + 7i2)/2r 



(85) 



and hence F^p ~ Hi. This immediately leads to 



P(ln.g) oc (4x 



2 

2'min 



2)e-T(r+")', u = ln{g/A). (86) 
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FIG. 22: Mean value (Ing) and variance var In g as a func- 
tion of 7i2 in the strongly insulating regime, calculated from 
analytical distribution 18311 . As expected, var \ng ~ —(Ing) 
for 712 ^ 1 but var Ing ~ — 2(lngi) for 712 — 1. Note that 
both mean and variance depends only weakly on 712 when 712 
is small (which is always the case when the system size L is 
large) . 



We do a saddle point analysis of Eq. to obtain 
(Ing) and var(ln(7) as a function of T. In order to il- 
lustrate the difference between QID and 3D insulators, 
we will keep the general expressions without using the 
condition 712 ^ 1. The free energy can be written as 



1 



approx 



{x - -)- - ln(2a;2,i 



-x); a;=:ln(l/.g). 

. 

The saddle point solution of the mean Xm = (In 5) is ob- 
tained from F'{xm) = where the prime denotes deriva- 
tive with respect to x. Denoting 



Xm = 1/r - y 
this gives y{y — 712/r) ~ 2/T, leading to 



y 



(211)2 

^ r ^ r r 



The variance Sx can be estimated from 1/F"(a;„ 

1 



6x 



4 



(88) 



(89) 



giving 



(90) 



In the limit 712 
tors, y \/2/T. 




-> appropriate for our 3D insula- 
On the other hand in the QID limit 
Thus compared to the QID result 
(Ing)Qi'D ~ i/r, the 3D result is shifted by Sim- 
ilarly, compared to the QID result var(lng) — 2/r, the 
3D insulators have a much sharper distribution, with half 
the variance l/F. Both results agree with numerical data. 
Although our model is not in general valid for 712 1 
because of our neglect of the interaction terms, in the in- 
sulating limit the interaction terms are negligible and it is 



FIG. 23: Skewness as a function of 712 in the strongly insu- 
lating regime. As expected, skewness is zero for 712 = 1 (QID 
limit). 
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FIG. 24: P(ln g) obtained from analytical formulas (Eqs. 
1831 and IMIl for F — 0.03 and various values of 712. 



useful to see how the mean and the variance changes with 
712 as it is changed from the 3D limiting value of zero 
to the QID limiting value of unity. Since 712 ~ C/2i, 
for a given disorder this will correspond to starting from 
a cubic sample of width L = ^ £, (where Lz is the 
length) and decreasing the width to L < ^ to reach the 
QID limit. 

Figure [221 shows how the mean and the variance 
changes with 712 according to Eqs. (|88|l and H90I) . 

It is not possible to obtain a simple formula for the 
skewness ((Ing— (ln5))'^)/[((ln,g— (Ing))^]'^/^ except that 
in the limit F — > it approaches a number of order unity. 
Direct evaluation of the quantity as a function of 712 is 
shown in Fig. (23 which shows that for a given disorder, 
the skewness starts from zero in the QID limit, as is 
well known, but saturates to a finite value (depending on 
disorder) in the 3D limit. It shows that the distribution is 
never log-normal for 3D insulators. It also shows that the 
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FIG. 25: Conductance distribution for 3D insulators obtained 
from direct numerical simulation for W = 16, L = 26 (cir- 
cles), and from Eq. ^ for T = 0.063 and 712 = T/2 (solid 
line). Both have the same mean value (ing) ~ —10.6. Dashed 
and dotted lines show Eq. I86II with a::2min = l/2r + 3/4 and 
2;2min = l/E, respectively, with F = 0.063. Shown are also 
numerical data for W = 26 and L = 10 (triangles). Similar 
agreement is obtained for other values of (Ing) if F is used 
as a free parameter; see Ref. [l5| for the case (Ing) ~ —12.6 
fitted with F = 0.054^. 



distribution P(ln g) is almost independent of 712 provided 
that both 7i2 and T are smaU. This explains why the 
distribution shown in Fig. 1251 does not depend on the 
ratio r/7. 

Finally, Fig. [21 shows how the entire distribution 
P(ln g) changes as a function of 712, from a sharp, skewed 
form for 712 ^ 1 to a broad Gaussian form for 712 = 1. 



B. Comparison with numerical data 

When comparing our theoretical prediction, Eq. (|S^ . 
with numerical data, we distinguish two cases: (1) for 3D 
systems, both F and 712 decreases ~ 1/L, while the ratio 
r/712 does not depend on L. (2) for systems x L^, 
r ^ l/Lz while 712 ^ l/L does not depend on Lz (see 
fig. |SJ). Therefore, contrary to 3D geometry, r/712 ^ 1 
for Lz ^ L. This different behavior of parameters F and 
712 explains difference between the shape of p(ln 5) in 3D 
and QID strongly disordered systems, shown in figures 
EUandEi 

Figure 1^31 shows Eq. compared with the results 
from direct integration of Eq. H83I) , both compared with 
numerical results based on Eq. (|11|) . For the analytic 
curves, we chose F = 0.063 and 712 = F/2 to have the 
same (Ing) as in the numerical case^S. Note that using 
the QID result 7^^^^ — 1 gives X2min ~ 1/r, leading to 
a log-normal distribution (see dotted line in Fig. H25|l ). 
As shown in [l5j |. variance and skewness calculated from 
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FIG. 26: Conductance distribution for insulating samples of 
the size x Lz {L = 10) In contrast to 3D system (fig. 
I25|l . parameter F/712 < 1. From numerical data, we have 
F = 0.3, 712 = 0.22 for W = 16, and F = 0.133 (0.102), 
712 = 0.32 (0.36) for systems with W = 14 and L^ = 30 (40), 
respectively. Solid lines show analytical result, Eq. 1831 1 with 
712 and (In g) as given form numerical data. 



direct integration of Eq. H83|) compares well with nu- 
merical results, consistent with saddle point results from 
Eq. PI). 

Figure compares theoretical formula, Eq. (|83|l with 
numerical data for p(hig) for insulating samples L^ x Lz- 
While F is small, decreasing as '-^ '^/Lz, 712 does not 
depend on Lz and is constant for fixed L. Consequently, 
ratio 712 /F ~ Lz/ L increases with increasing Lz- 

Both figures|211and[5B|show qualitative agreement with 
numerical data and theoretical model. Quantitative dif- 
ferences between Eq. (|83f) and numerical results have ori- 
gin in our simplified model Eq. (|19|l . which still overesti- 
mates the strength of the interaction for higher channels. 

It is important to note that in both figures 1251 and 
1261 we compared numerical data with theoretical model 
with the same mean (In g) . This is consistent with scaling 
theory of localization since there is only one parameter - 
for instance (In (7) - which determines p{\ng) completely. 
In order to make sure that the analytical model has the 
same (In 17) as the numerical data for a given disorder, we 
used F as a free fitting parameter. Of course we could use 
Eq. (20) to obtain F independently for a given disorder. 
However, in order to do that we will need a good estimate 
of the mean free path I. This is difficult in the strongly 
disordered regime because the mean free path defined as 
the decay length of the single particle Green's function is 
actually smaller than the lattice spacing in the strongly 
disordered regime*^^, and our numerical model does not 
allow us to obtain such small lengths with good accuracy. 
While independent calculations of the mean free path are 
available for cubic systems below critical disorder— (e.g. 
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FIG. 27: Mean value (Ing) as a function of l/F. Open sym- 
bols: analytical result, Eq. 1831 . Data confirm linear de- 
pendence (Ing) oc 1/r for small F (strong disorder). For 
smaller disorder (larger F) analytical model is less accurate. 
Data confirm that (In g) does not depend on 712 provided that 
both F and 712 are small. Full symbols: numerical data with 
F~^ = LKii/£. We estimate F for a given disorder using 
the mean free path £ obtained from data in fig. 19a {£{W) 
is shown in inset) and the relation K\\L = xil. Numerical 
(Ing) is an unambiguous function of F, as required in the an- 
alytical model. Small deviations are due to finite size effects 
which affect actual values of parameters xi and Ku. 



i — 0.234 for = 15 in the isotropic case), there is no 
data available in the insulating regime. We therefor use F 
as a free parameter. Nevertheless, as a consistency check, 
we estimate I from a plot of K\\L vs x\ (fig. 19a), where 
the slope should give the mean free path (equivalently, 
we could identify F"^ with numerical value of x\). The 
results are plotted as the inset in fig. |27|to show the mean 
free path as a function of disorder. Using this result, we 
can estimate the value of F corresponding to the disorder 
= 16 used in fig. 25. We find that liy^ = 16) « 0.48 
and F = 0.070, which is close to the fitting value 0.063. 
This shows that while we can not obtain F accurately 
enough in our present numerical scheme, the fitted values 
are consistent with our crude estimates. As a further 
consistency check, we use the above estimate of the mean 
free path to plot in fig. 27 the F dependence of (Ing). It 
shows that (In g) is indeed an unambiguous function of 
F, as required by the theory. 

Finally, we have checked the effects of fluctuations of 
fell on P{\Mg) by integrating the conductance distribu- 
tion in fig. 25 over the distribution P(fcii) (fig. 2J. We 
find that the effects are negligible. 



VIII. BEYOND THE INSULATING LIMIT 

From numerical simulations^i^ we know that the criti- 
cal regime in 3D is also dominated by only a few eigen- 
values Xi 3> 1. However, since 712 is neither nor 2, it 




FIG. 28: (a) Probability distribution P(ln g) for F ~ 1 where 
critical regime is expected. The distribution agrees qualita- 
tively with numerical data for the critical regime. Figure (b) 
shows detail of the distribution shown for F = 0.35, 712 = F/2 
and for F = 0.6, 712 = 1. It shows that there is indeed a non- 
analyticity in the distribution close to In g = with position of 
the non-analyticity at Injr > 0, in agreement with analytical 
results-. 



seems that we may not be able to use the free fermion 
(77 = 0) solution of Eq. H24|) to obtain the distribution 
of the transmission levels. However, as shown in |32j . 
the solution is independent of the strength of the inter- 
action 77 in the strong disorder regime characterized by 
3> 1. This means that our solutions might be used, al- 
beit only qualitatively, even near the critical regime. We 
show in Fig. (|28|) the distribution P(ln 5) for F ~ 1 which 
is expected to be near the critical regime. It agrees qual- 
itatively well with numerical results at the critical point, 
including a discontinuity in the slope near g = \. It is 
known from analysis of the QID systems^ that separating 
out an additional level helps to study the non-analyticity 
near g — \. We therefore expect to obtain better results 
near the critical regime by separating out an additional 
level. 

Finally, we show in Fig. 1551 the distribution P{g) for 
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FIG. 29: Analytic result for the metallic regime (F > 1 and 
7i2 = 1). Data confirm that the distribution P{g) is Gaussian 
with variance var g — 0.064 for F — 1.5 and var g — 0.09 for 
F = 3.5 which is comparable to the UCF value. 



r > 1 and for 712 = 1 which corresponds to the metal- 
he regime. Although we do not expect that our ap- 
proximate formula works quantitatively for the metallic 
regime, Eq. H83|l gives, for this choice of the parame- 
ters, a Gaussian distribution of the conductance. Also 
the width of the distribution qualitatively agrees with 
the universal conductance fluctuations^ in this regime. 
This shows that our simple model already captures the 
essential qualitative features at all strengths of disorder. 



IX. SUMMARY AND CONCLUSION 

We systematically analyzed the length and disorder 
dependence of the matrix K to check if the generalized 
DMPK equation proposed in ^T^l is valid in three dimen- 
sional systems at all strengths of disorder. 

We studied the matrix K in detail. The goal was to test 
the assumptions on which the generalized DMPK equa- 
tion was derived and to construct a simple analytically 
tractable model for K which captures all the important 
qualitative features. In particular, since QID systems 
have been studied in great detail, we looked for any ma- 
jor qualitative differences in the structure of K between 
QID and 3D systems. 

We find that to a good approximation, the generalized 
DMPK equation remains qualitatively valid for any disor- 
der. We also conclude that to a good approximation, we 
can use only two parameters, Kaa ~ Kn and ^ab ~ 712, 
to characterize the qualitative changes in transport at dif- 
ferent strengths of disorder in different dimensions. We 
find that although fluctuations in fcn at strong disorder 
are large (non self-averaging), the effect of these fluctu- 
ations on P(ln g) is negligible. More importantly, we do 
not have an independent way to estimate the mean free 
path to obtain F = i/ (LzKn), but all qualitative fea- 




Ing 



FIG. 30: Typical form of the distribution p(lng) in 2D disor- 
dered systems compared with Gaussian distribution with the 
same mean value and variance. (In g) — —20.6, var Ing = 23.5 
and skewness is 0.251. Figure indicates that the deviation 
from the Gaussian distribution exists in 2D as well. 

tures of the entire distribution P{\ng) is obtained cor- 
rectly once an effective F is used as a free parameter. We 
also find how these parameters depend on disorder and 
show their unambiguous dependence on the localization 
length. This is important since it indicates that the intro- 
duction of new parameters does not necessarily invalidate 
the single parameter scaling theory of localization. 

We have also shown that the matrix K contains infor- 
mation about the Anderson transition. The scaling of 
the parameters LKn or 712 clearly identifies the critical 
point which agrees with numerical results. 

We then concentrate on the strong disorder limit where 
our numerical results allowed us to construct a sim- 
ple one-parameter model of the matrix K, containing 
r = e/KnL, and 712 - 2X^2/ K^, with r/712 = 2. 
By varying 712, we show how one can go from a QID 
^ L) to a truly 3D <^ L) system in the insulat- 
ing regime, which clearly shows the difference between a 
QID and a 3D insulator. We then use the model to ob- 
tain the full distribution P{g) which agrees qualitatively 
with numerical results. 

It is indeed remarkable that even though the gener- 
alized DMPK equation H22|) neglects fluctuations in kab 
and the model Eq. (19) neglects the index dependence of 
Kab, the theory still captures all the essential features of 
length, disorder as well as dimensionality dependence of 
the entire conductance distribution and provides in par- 
ticular a simple understanding of the 3D distribution at 
strong disorder, which is qualitatively different from a 
log-normal distribution in QID. At the same time, our 
numerical studies suggest that having an independent es- 
timate of the mean free path could provide a more quan- 
titative description of the conductance distribution in 3D 
at all disorder. 

We emphasize that there are large differences between 
QID and higher dimensions. In QID defined in PIsll^fTsl. 
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disorder is always weak enough to assure that the local- 
ization length ^ 3> i where L is the transverse dimension. 
The QID insulator corresponds to the weakly disordered 
systems of length ^ C It is this length-induced insu- 
lating behavior that is described by the DMPK equation. 
This is different from localization in 3D which occurs at 
strong disorder, where ^ is much less than both and 
L. This difference is clearly reflected in the matrix K, 
where Table I summarizes how the scale dependence of 
Kii and 712 depend on disorder in 3D. In contrast, K 
in QID is independent of disorder. Our model recovers 
all the peculiarities of the 3D localized regime: we found 
that the distribution P(ln g) is narrower than in QID and 
possesses non-zero skewness. 

Although we concentrated on 3D systems, general con- 
siderations about the properties of the matrix K in the 
insulating regime should be valid in any dimension d > 2. 
In particular, the assumption that Kn is a nonzero L- 
independent constant in the localized regime is correct 
independent of dimensionality. Therefore, we expect that 
our theory of the insulating regime is valid for any d > 2. 
Consequently, the distribution P(ln g) is not Gaussian in 
the localized regime in any d > 2, although the nature of 
the deviation from the Gaussian form might be dimen- 
sionality dependent. In fact we expect P{\ng) to be dif- 
ferent from Gaussian distribution even in 2D^^. This ex- 



pectation is supported by Fig.l30lwhich shows P{\ng) for 
a 2D square system obtained numerically using Eq. ITTI 
As discussed in the paper, deviations from the Gaussian 
form are due to the changes of the spectrum of parame- 
ters X. Although the contribution of the first channel to 
the conductance is dominant, higher channels do influ- 
ence the statistical properties of the smallest parameter, 
xi. This effect was probably not considered in previous 
analytical works which predict Gaussian distributions of 
In g in the insulating regimes in dimensions d = 2 + e~. 

When applied to the critical regime, our theory recov- 
ers typical properties of the conductance distribution, 
including the non-analyticity of the distribution in the 
vicinity oi g = 1. Although our present results are only 
qualitatively correct in the critical regime, we believe that 
the method developed in the paper represents a good 
starting point for further development of the theory of 
Anderson transition. 
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